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Abstract. We discuss the prospect of extracting the values of the mixing parameters S and 613 
through the detection of cosmic neutrinos in the planned and forthcoming neutrino telescopes. 
We take the ratio of the //-track to shower-like events, R, as the realistic quantity that can 
be measured in the neutrino telescopes. We take into account several sources of uncertainties 
that enter the analysis. We then examine to what extent the deviation of the initial flavor 
composition from w e : : w T — 1 : 2 : can be tested. 



The neutrino mixing parameters can be extracted from cosmic neutrinos data according to 
the following argument [lj: assume that the flavor ratio of neutrinos at the source is w e : : w T ; 
after that neutrinos travel the large distances between the astrophysical sources and the Earth 
the oscillatory terms in the flavor transition probabilities average out such that the flavor ratio 
at the detector will become 

w a P ae : w a P afM : w a P ar , (1) 

a. a. a. 

where 

P af 3 = P{V* -» Vp) = P(v a -^Up)=Y, \Uai\ 2 \U P i\\ (2) 

i 

and U a i are the elements of the neutrino mixing matrix. In a wide range of models the flavor 
ratios at the source are predicted to be w e : : w T = 1 : 2 : 0. Thus, by measuring the 
flavor ratio at Earth, one can derive the absolute values of the mixing matrix elements which in 
principle yield information on the yet-unknown neutrino parameters #13 and 5. 

IceCube [2] and its Mediterranean counterparts ( such as KM3NET [3]) can basically 
distinguish only two types of events: 1) shower-like events; 2) //-track events. It is possible 
to derive information on the flavor composition of neutrinos by studying the ratio R =number 
of //-track events/number of shower-like events [3]. In the analysis in this paper we consider 
neutrinos with energies 100 GeV < E v < 100 TeV, where the upper and lower limits come from 
the absorption of neutrinos in Earth and the energy threshold of detection in neutrino telescopes, 
respectively. 

Two sources contribute to the //-track events: (i) Charged Current (CC) interaction of Uu 
or producing // or //; (ii) CC interaction of v T and 9 T producing r or f and the subsequent 
decay of r and f into // and //. In the literature, the contribution of v T (via v T — > r — > //) to 



/Li-track events has been overlooked but to study the effect of #13, one should take into account 
such sub-dominant effects (see the Appendix of [5] for details). Three types of events appear as 
shower: i) the Neutral Current (NC) interactions of all kinds of neutrinos; ii) the CC interactions 
of v e and v e ; iii) the CC interactions of v T (i7 T ) and the subsequent hadronic decay of r (r). 
Details of the event rate calculation for each case in both //-track and shower-like events can be 
found in [5]. 
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Figure 1. (a) The dependence of R on sin 2 #13 for different values of the spectral index, a. The thicker 
lines correspond to 8 = tt and the thinner ones correspond to 5 = 0. In drawing this figure we have set 
NvJN Vc = 0.5 and (A/^ + M Vfl ) / (A/p e +M Ve ) = 2. The input for 6 12 and 823 are set equal to the best fit 
in [6]. The vertical line at 0.041 shows the present bound at 3a [6] (b) Points in the (a, S) space consistent 
with R = 2.53 ± AR. True values of the (a, S) pair are (2, tt/2). Points displayed by dots, triangles and 
crosses respectively correspond to AR/R = 7%, AR/R = 1.5% and AR/R = 1%. To draw this figure we 
have varied sin 2 13 G (0.028, 0.032), sin 2 9 12 G (0.30, 0.34), sin 2 6 23 G (0.47, 0.53) and N„JN Vc € (0, 1). 

Several input parameters enter the calculation of R and their uncertainties induce imprecision 
in the calculated value of R. We consider the following uncertainties in the calculation of R: i) 
for the energy spectrum of neutrinos we assume a power-law spectrum dF U0 /dE U/3 = N Uf3 E~®. 
Neutrino production through the Fermi acceleration mechanism of particles in the source predict 
the value of the spectral index equal to a = 2, but taking into account non-linear effects in 
the acceleration mechanism results in spectral index values a € (1,3). It is shown in [1] by 
assuming the flux E^,dF v / dE u = 0.25 GeV cm -2 sr _1 yr _1 , after one year of data-taking a can 
be determined with 10% uncertainty. For the normalization factor J\f u „ it can be shown that: 
= H Vlk and N„JN Ve G (0, 1). 

Fig- d-a) shows R versus sin 2 #13 for cos 5 = ±1 and various values of a. As seen from the 
figure when 5 = 0, the sensitivity of R to s 2 3 is very mild and less than 2 %. That is while for 
cos 5 = —1, the sensitivity to s 2 3 is about 10 %. As seen from the figure, even for cos<5 = — 1, 
the sensitivity to s 2 3 can be obscured by the 10 % uncertainty in a. However, for s 2 3 > 0.02, 
the bands between a = 2.2 and 1.8 for cos 5 = 1 and cos<5 = — 1 have no overlap. This means 
that for s 2 3 > 0.02, 10 % precision in a is enough to distinguish cos 5 = 1 from cos 5 = —1. 

Fig. {H-b) addresses the question that whether it will be possible to extract the value of 5. 
Drawing the plot, we have assumed that R will be found to have a typical value of 2.53 with 
an uncertainty of AR/R. This value of R can be obtained by taking maximal CP-violation 
(5 = tt/2), sin 2 (9i 3 = 0.03, w e : : w T = 1 : 2 : 0, a = 2 and MpJM Ve = 0.5 . We have 
looked for solutions in the 5 — a plane for which R = 2.53(1 ± AR/R), varying the rest of the 
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Figure 2. Points in the (w e ,w T ) plane consistent with R = R±AR. The ratios are normalized such that 
w,, = 1. The true values of {w e , w T ) are denoted by -jc. Points displayed by dots and triangles respectively 
correspond to AR/R — 7% and AR/R = 1%. In drawing this figure we have varied sin 2 #13 £ (0, 0.003), 
S G (0,2tt), a S (1.8,2.2), and AfpJAf^ G (0,1). Drawing Fig. (a), we have taken R = 2.53 which 
corresponds to the standard picture with w e = 1/2 and uv = 0. In case Fig. (b), we have set R = 3.2 
which corresponds to the stopped muon scenario with w e = w T = 0. 



relevant parameters in the ranges indicated in the caption of Fig. (JTJ-b). The regions covered 
with dots, little triangles and crosses respectively correspond to 7%, 1.5% and 1% precision in 
the measurement of R. As seen from the figure, with AR/R = 7%, 5 cannot be constrained. 
In fact, any point between the vertical lines can be a solution. The figure shows that reducing 
AR/R to 1% (but keeping the rest of the uncertainties as before), some parts of the solutions 
can be excluded. In particular, the region around 5 = tt will not be a solution anymore. Notice 
that along with 5 = ir/2, 5 = is also a solution. This means that despite maximal CP-violation, 
the CP-violation cannot still be established. 

Taking into account the relevant uncertainties in the input parameters, we look for values of 
w e : Wn '■ w T that are consistent with R = R±AR. To perform this analysis, we take #13 = 0. In 
Fig. ([2]), we consider two possibilities: (i) the standard case with w e : : w T = 0.5 : 1 : leading 
to R = 2.53 (see Fig. [2]-a); (ii) the case of stopped muons with w e : : w T = : 1 : yielding 
R = 3.20 (see Fig. [2]-b). From these figures we observe that with a precision of AR/R = 7%, 
these two scenarios can be easily discriminated. These two can also be discriminated from the 
scenario in which the neutrino production mechanism is n — > pev e {i.e., w e : : w T = 1 : : 0). 
When we restrict the analysis to w T = {i.e., the case without exotic neutrino properties) from 
these figures we observe that the measurement of R stringently constrains w e : w T which in turn 
sheds light on the production mechanism. However, once the assumption of w T is relaxed, a wide 
range of w e : : w T can be a solution. For example, the exotic case of w e : : w T = : : 1 
leads to the same value of R as the stopped muon scenario. 
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